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Discussion 


Several  recent  experiments  have  shown  features  that  indicate  the  presence 
of  an  Ha"  ion,  Schultz1  proposed  a  state  of  Ha"  in  the  formation  of  E  and  H" 
lay  electron  impact  on  hydrogen  gas.  Khvostenko  and  DtJiel^shii23  proposed  a 
similar  mechanism  in  electron  impact  measi*rements  on  Es,   i.e.,  H2  +  e  — >  (E2""5 
H  +  E~.  They  also  detected  a  negative  ion  of  mass  two  by  mass  spectrometer213 
upon  "bombarding  a  mixture  of  water  and  antimony  gas  with  electrons.  In  ref. 
(2b)    is  explained  why  this  was  felt  to  be  Ea~  with  a  lifetime  long  enough 
to  be  observed,  rather  than  the  other  possibility,  D".  In  astrophysics3  it 
has  been  suggested  that  Hs~  occurs  in  the  solar  corona  and  dissociates  into  E 
and  E~\ 

The  latest  of  these  have  been  some  electron  scattering  experiments  by 
Simpson  e.ndi  Kuyatt,4*5  where  a  highly  monochromatic  electron  beam  was  seatterec. 
off  hydrogen  gas.  A  plot  of  current  vs.  energy  was  made,  and  resonances  correE- 
ponding  to  a  set  of  vibrational  levels  i:ere  observed,,  From  these  levels  the 
spectroscopic  constants  coe,  a>exe,  and  the  approximate  E0  were  calculated.  The 

spectroscopic  constants  are  very  similar  to  those  of  E2,  ^If^  or  ^j, -however  the 
threshold  energy  was  lower  than  both,  thereby  suggesting  the  possibility  cf  a 
resonant  state  of  Eg"  consisting  roughly  of  an  inner  a   electron  and  two  outer 
electrons,  with  symmetry  2<rg«*  (Note  that  the  outermost,  planetary  electron  would 
not  greatly  affect  the  field  in  which  the  nuclei  move.) 

The  predominant  common  feature  of  all  these  experiments  is  that  they 
are  of  a  scattering  nature,  and  thus  the  observations  could  be  of  a  real  bound 
excited  state,  or  of  a  resonance  in  the  continuum. 

Previous  papers  concerning  theoretical  calculations6  on  Eg"  show  that 
there  is  no  bound  state  with  an  energy  lower  than  that  of  the  Eg  X  1Zg  state 
of  the  hydrogen  molecule  with  an  electron  at  infinity.  In  the  computations 
of  (6)  many  configurations  were  tried,  with  complete  variation  of  non-linear 
and  linear  parameters,  and  in  all  cases  the  variation  method  itself  prod\<ced 
the  equivalent  of  a  hydrogen  molecule  with  an  electron  in  the  continuum. 

Thus  it  was  decided  to  look  for  a  bound  state  of  a  different  nature,  i.e., 
following  the  suggestions  of  experiment,  a  state  with  one  inner  c   electron  and 
two  outer  y  electrons,  such  that  essentially  the  wave  function  would  resemble 
the  trapping  of  an  electron  in  a  tr  orbital  by  a  1T^l  or  3%x   hydrogen  molecule. 
This  appears  physically  resonable  since  an  incoming  electron  would  "see'  a  a 
electron  aboutAwo  nuclei,  and  a  Itr  electron  perpendicular  to  it,  or  therefore, 


*    L*  and  E"  are  essentially  degenerate  for  this  discussion .  This  can  be  seen 
from  the  small  overlap  of  the  (S,r\)   parts  of  lirand  3?r» 


a  net  attractive  charge  whicl:  was  greater  than  zero,   since  the  excite-     o- 
gonal  outer  electron  would  not  screen  for  the  incoming  one-  Also,  in  br 
the  two  electron  (lo  lur)   and  the  three  electron  (lor  ."  ,r  ^w)   -systems,  the  electn 
are  highly  un cor related.  The  ifciportance  of  this  will  be  di.  seussec,  :  3 

Note  the  difference  between  this  and  electron  capture  by  a  hydro* 
atom  to  form  H~o  In  atoms  aii  incoming  electron  can  i?iduee  polarization 
the  hydrogen  molecule,  most  of  the  electron  density  must  remain  between 
nuclei  to  form  the  bond,  and  therefore  cannot  be  polarized  nearly        as 
two  separate  atoms  by  the  incoming  electron.  This  gives  another  phys 
of  the  reason  a  ntate  like  the  1T,g  state  of  hydrogen  cannot  capture  a 
but  states  like  LTTU  and  ^H^  night  be  likely  to. 

Use  of  a  trial  variation  function  of  this  type  could  lead  bo  two 
results.  First,  but  highly  unlikely  due  to  experimental  e-vidence  tad  p:   f 
theoretical  calculations,  there  could  be  a  true  bound  state  which  woo     a 
solution  to  the  Schroedinger  aquation,  with  a  true  variation  princ.lpl 
energy  less  than  that  of  the  ground  state  of  the  hydrogen  volecule.  Or 
the  lc  ly  configuration  could,  cause  a  narrow,  but  very  deej:  potential 
the  incoming  electron,  such  that  there  is  a  resonance,  comr.osed  of  a  h 
localized,  relatively  long-lived  combination  of  continuum  states  (wave  packe 

The  problem  then  is  to  construct  such  a  highly  localized  wave  ] 
from  continuum  functions  having  a  lifetime  long  enough  to  ze   observable  in 
other  words,  with  a  low  probability  of  a  radiationless  transition  fron  ,he 
resonant  state  to  the  hydrogen  molecule  plus  an  electron  of  proper  energy  in 
the  continuum. 

Using  the  physical  arguments  above,  a  configuration  which  will  hopefully 
have  a  deep  potential  well  for  the  scattered  election  and  a  smell  ove:.    rith 
the  Ra  1F.g  +  e  irust  be  constructed.  A  preliminary  calculation  was  made  'ith  a 
trial  function  with  configuration  loirlv  .  Even   this  rough  calculate 
the  suspicion  teat  there  is  no  bound  discrete  state  with  an  energy  less  ;;hsn  that 
of  H2  1rg  +  e,  The  second  possibility  then  becomes  the  one  which  sua     In- 
vestigated. Certain  difficulties  arise,  for  there  it;  no  longer  s  varl 
theorem  and  no  rigorous  minlaum  principle.  How  does  one  select  the  pe 
of  a  wave  packet  and  calculate  the  energy  under  thes.e  circimsstances? 
raents  will  be  given  to  show  the  validity  of  an  approximate  variation  princi 
They  are  based  en  the  fact  that  the  syBtezr  is  to  a  large  d-sgree  uneorrff  Inted, 
with  weak  interactions  between  electrons,  and  that  the  des->  potential  causes  ? 
resonance  packet,  which  is  a  function  of  time,  and  this  ti?ae  dependence  may,  to 


a  first  appro 

Heither  of  these  argument    i  be  prove?;.  ie   subsequen 

lion  and  results  show      x>th  to  be  physically  quite  valid. 

The  variation  j        requires      all  bound  states  he  lal 

to  all  lower  energy  states  >f  the  same  syiBRetry.  In  the  j:  resent  sa     he 
packet  ap        ng  a  bound  state  must  be  orthogonal  to  all  zZg   3:  :    of 
Hg  1 Eg  plus  an  electron  at  Infinity.  Al      rigo; :     bhogonalitj 
possible  since  both  s"       a  2Zg,  effective  orthogon     [zero  01 

ible  if  the  packet,  is  constructed  so     i.wo  inner  electrons 

iogonal  to  He  Ion  is  restrl  e  tve 

inner      ras  must  remain  la,  liru.  i.  efore 

orthogonal  to  H,»  1Zg.  Orthogonality  is  also  expected  :. 
electron  in  the  lover  state  overlaps  negligibly  with  the  localiz 
packet 0  This  has  been  verified  by  direct  calculation  of  >er 

of  non-equal  energy  cases.  If  the  pack  re  to  overlap  with  1 
state  it  would  have  to  be  those* of  equal  energy.  The  foj  owing 
presented  to  account  for  this  case.  It  has  been  not: 

good  agreement  with  experiment,  the  in-     Lvity  of  the  art  :~ 

trcn's  parameters  to  the     's  parame,  'ting 

molecular  orbitals  indicate  that  the  electron  correlation  is  vc         If 
this  is  the  case,  l/ria  can  be  consider;     a  perturbation  and  a  zero 

that  the  overlap  is  zero  can  be  given.  Higher  o:fder  perturbation  theory 
will  give  negligible  overlap  if  there  ii    avergence"  by  zero   order. 
"^corrJ^  *s  3MaH  then  the  configuration  used  is  a  product  of  three  llzr   functions. 

packet  i3  lag  liru  x  3hfu_1J  *fee  cont     is  la  lw, 

rly  orthogonal  by  the  orthogonality  of  Ha+  orbitals.  There! 
bat  ion  theory,  weak  correlation  implies  that  a  res:       ariatioi 
to  different  H^r  orbitals)  allows  a  true  atinimum  prin-.     ?ith  which 
late  the  "best"  energy,,  To  verify  this  it  is  noted  in  Holbein,  anc 

ork  on  helium7  very  reasonable  results  are  obtained  when  th« 
Ls  applied  to  such  vncorrelated  packet  states  a 
(2p3p)  xPg,  (2p3d)  1DU,  anc  (2p3d)  9Da> 

(2»)*  1S  state7**  in  H"  has  a  lower  energy  than  that  o  state 

of  H".  In  helium,  the  reverse^  expected  result  is  found,,  It  is  felt  hy   the 
authors  that  the  result  on  H"  is  in  error  because  in. the  (i 
tion  is  not  small  and  no  approximate  variation  principle  isolds . 


"being  lower  than  (Is)2  is  unphysical  and  contrary  to  the  more  extensive  calcu- 
lation of  Pekeris8a  where  he  shows  only  one  bound  state  exists.  The  (2s)2  1S 
problem  in  helium  is  also  not  completely  settled. Qa'fe  This  is  felt  to  be:  due 
to  the  correlated  nature  of  the  system  and  the  consequent  failure  of  the 
variation  method. 

Another  justification  of  the  variation  principle  is  that  since  this  is 
a  resonance,  the  energy  calculated  by  variation  is  not  a  discrete  eigenvalue, 
but  the  average  of  a  spread  of  energies.  This  average  energy  must  be  a  function 
of  time,  since  the  packet  must  eventually  spread.  Now  there  is  no  explicit 
time  dependence  in  the  variation  function,  but  the  variable  parameters  regulate 
the  spread  of  the  wave  packet,  and  thus  the  parameters  themselves  may  "be  con- 
sidered as  time  dependent  and  tending  to  cause  the  resonant  wave  packet  to 
delocallze  as  time  increases. 

To  achieve  the  longest  lived  resonant  state  in  'the  restricted  region 
of  Hilbert  space  selected,  therefore,  by  the  uncertainty  principle,  &&&v  must 
be  minimized. 

This  can  be  stated  as  E^y  *  %&■?(&!),   where  a±  are  the  variable  parameters, 
and  qjjl  *  «i(t).  Then  to  minimize  AE, 

^w  .  o  -  S^t^  (I.l) 


5t   -    u  -  $&£ 

But  the  variation  method  causes  dEaV/<&i  =»  0  for  all  o^  therefore  min- 
imizing AEgv  and  maximizing  At.  This  in  essence  selects  the  best  (longest 
lifetime)  wave  packet  possible  in  the  region  of  Hilbert  space  considered,,  Note 
that  the  variation  alone  will  not  give  a  long  lifetime  state,  but  will  only 
give  the  packet  of  longest  lifetime  among  the  class  chosen. 

The  use  of  the  variation  principle  has  now  been  physically  justified. 
Thus  the  question  becomes:  is  the  state  which  has  been  found  of  long  enough 
lifetime  to  be  observed?  An  appeal  to  the  ultimate  agreement  with  experiment 
could  be  made,  with  the  claim  that  the  state  clearly  exists  long  enough  to 
vibrate.  However,  it  is  preferable  to  give  a  physical  argument  based  on  the 
weakness  of  the  correlations  and  the  resulting  perturbation  method  to  indicate 
a  long  lifetime o 

To  demonstrate  this,  consider  the  standard  auto-ionization  picture  of 
bound  states  in  the  continuum.9  The  argument  will  be  outlined  here  for  helium 
because  it  is  easier;  it  will  be  clearly  extendable  to  the  H2~  case.  If 


iA'ia  ■  Vcorr,  is  small  the  zero  order  wave  functions  are  products  of  two 
hydrogenic  orbitals.  The  functions  and  levels  fall  into  three  classes  with 
overlapping  energies: 

(a)  E  a a g    both  one-electron  functions  bound 

1    k©*^ 

(b)  B  * 1  +  -5-    one  one -electron  function  bound  and  one  continuum 

111     2 

kia   kofi 

(c)  E  «   L£-  +  -*-    both  continuum 

When  nx   and  n2  are  large  in  (a)  and  nj.  small  in  (b),  a  k  can  be  found  such 
that  Ecorit0  m  Ebd.>  i«e«>  there  exists  a  bound  state  in  the  continuum,  A  non- 
radiative  transition  from  the  bound  to  continuum  state  of  equal  energy  is  given 
by  the  probability  of  auto-ionization  (see  ref.  9,   and  ref.  1&),     Therefore  the 
transition  probability,  w,  is  proportional  to 

[<Wid  j  jLj  J  cont.)J  2  bOfa.  -  EcontJ 

The  lifetime  of  the  bound  state  is  proportional  to  w"1. 

Now  if  the  correlation  is  small,  an  assumption  upon  which  the  equation 
for  w  was  derived,  the  first  order  correction  to  the  bound  state  function  due  to 
correlation  is  given  by 

(Ho  -  Bbd.)  *i  +  (^  -  Ex)  fto±9   -  0 

If  ti  is  expanded  as  usual  on  all  other  states  of  the  system,  and  substi- 
tuted back  into  the  perturbation  equation,  then  multiplied  by  tycont,  an<i  integrated,, 
the  result  is,* 

^boundj  ~-|  cont.)  +  Ofad.  ^  (%,  -  Econto)  »  0 
ia  cont. 

Now  when  E^  ■  Econt.  the  matrix  element  is  to  first  order  zero  and  w  is 
then  of  the  form  x8(x),  x  «  B^  -   Econt.'  an<i  therefore  approximately  zero. 
Higher  order  perturbation  theory  modifies  this,  but  if  the  perturbation  is  really 
small,  in  a  negligible  way. 

In  support  of  this  argument  it  is  pointed  out  that  the  types  of  resonant 
states  that  have  been  observed  are  generally  highly  uncorrelated  and  ones  in 
which  the  "resonant  trapped"  electron  sees  a  deep  well  due  to  the  rest  of  the 
systemo  It  is    exactly  this  type  of  state  for  which  successful  variation 
calculations  have  been  carried  out. 


*  In  this  expression,  the  continuum  normalization  factors  have  been  cancelled 
from  both  sides. 


II.  Wave  Function  and  Calculation  Methods 


Using  the  physical  arguments  above,  it  was  decided  to  use  a  one  con- 
figuration three-electron  variation  function,  having  ^g  symmetry,  such  that 
loss  of  the  outer  electron  would  not  he  ahle  to  yield  a  2Eg  state,  hut  only 
a  ^ITy  state.  The  elliptic  molecular  orbitals  and  calculation  methods  of  Earris 


(II. 2) 


and  Taylor  were  used.10 

For  one  configuration  the  variation  function  "becomes: 

V  W>  *  6*  t  E.  IL(P)  $(  Pr)   rft.  C5)  (II.1) 

where  Uj<i(P)  is  the  Jd  element  of  the  matrix  representation  of  the  permutation 
P  on  the  spin  eigenfunction  bases.  The  €^   are:11 

0^(5)-    AT4,   (a£A-|*a) 

The  spatial  part  is  the  product, 

XCT^*TT  0;^-)   +TT   &'<*<>  (H-3) 

where  the  0j  are  of  the  form, 

and  the  0^'  are, 

(-0      0;  (£,-*;  w^-y)  (II. 5) 

|,  tj  and  0  being  the  "standard  two-center  elliptical  coordinates.  The  b±f  <x± 
are  variable  within  the  computer  program  used,  but  n,m,v  are  input  parameters 
and  therefore  fixed  throughout  the  problem.  The  configuration  chosen  was; 

It  is  convenient  to  denote  these  in  terms  of  the  v±,   since  none  of  the  n,m,v 
change  and  all  m^  «  0.  Thus, 

X<?>=  ?U,(jJ_  +*'#&'  (n-7) 

Since  Em^  «  0  and  Zlv^l  ■  2,  there  is  no  sign  change  between  the  products, 
and  gerade  symmetry  is  ensured.  Also,  Sv^  «  0,  so  the  three  electron  product 
is  a  £  state. 

Note  that  if  the  outer  electron  is  removed  (placed  in  the  continuum), 

is  obtained,  which  is  aTT  state  (Zv  ^  0)  with  ungerade  symmetry  since  the  01/' 


7 

introduces  a  -  sign  between  the  products.  Thus  it  has  been  spatially  ensured 
that  the  1  configuration  variation  function  cannot  vary  itself  to  Ha  1Eg  +  e 
by  varying  6-1  to  a  very  small  number  which  has  the  effect  of  ejecting  the 
outer  electron  into  the  continuum. 

Looking  at  the  spin  portion,  it  is  also  noted  that  &yat.      can  "lose"  an 
electron  to  become, 

&Z»*^  s  £    tW-^^5)  (n.9) 

which  is  the  spin  eigenfunction  for  a  two  electron  singlet. 

The  one -configuration  variation  function  is  now  constrained  in  such  a 
manner  that  it  cannot  eject  the  outer  electron  and  become  anything  but  H2  ^n^ 
•5-  e,  and  the  resonant  wave  packet  is  suitably  localized  by  these  constraints. 


8 

III*  Calculations  and  Results 

In  order  to  obtain  a  potential  curve  and  the 
spectroscopic  constants  for  comparison  with  those  reported  by  experiment, 4  it 
^as  planned  to  calculate  five  energy  values  ab  nuclear  radii  corresponding  to 
scaled  roots  of  the  fifth  order  Chebyschev  polynomial  Ts(x),  where  x  ■ 
R  -  Ke/.3k026,   which  gives  separations  of  ±   .2000  and  ±  ,5236  au,  about  Rq, 
and  fit  these  to  Chebyschev  polynomials  To  through  T4.ia  The  spectroscopic 
constants  could  then  be  calculated  by  the  aethod  of  Dunham,,13 

Preliminary  calculations  with  a  trial  function  of  the  type  ltf(ttr)fi 
yielded  an  equilibrium  radius  of  1,9196  au.,  and  the  final  full  variation  of  the 
5  and  a  parameters  was  done  about  this  Re.  Upon  fitting  to  the  Chebyschev 
polynomials  and  minimizing  E(R)  with  respect  to  R,  Rq  was  found  to  be  1.9^9  au. 
instead  of  1,9196,  so  Dunham's  method,  based  on  a  polynomial  in  (R  -  ReV'Rey 
becomes  excessively  tedious.  Therefore,  the  polynomial  in  x  was  converted 
to  a  polynomial  in  R  and  the  vibration-rotation  constants  calculated  by  the 
method  used  by  Fraga  and  Ransil  in  ref .  14. 

To  note  the  effect  of  greater  interatomic  distances  three  other  R  values 
were  used.  Table  I  gives  the  potential  curve  and  parameters  of  the  wave 
function. 

Table  II  is  a  comparison  of  the  calculated  spectroscopic  constants  with 
those  reported  in  the  scattering  experiment. 

It  is  also  of  interest  to  note  in  Table  III  that  as  R  increases,  the 
potential  curve  approaches  the  curve  of  the  ^^  state  of  H2.15  This  can  be 
explained  by  inspection  of  the  parameters  of  the  three  electrons,  At  R  «  10,0  au, 
there  is  essentially  a  Is  hydrogen  atom  and  a  Ip  hydrogen  atom  widely  separated, 
with  an  electron  in  between.  In  H2  11J,U  at  this  internuclear  distance,  the 
expected  configuration  is  a  pair  of  hydrogen  atoms,  one  in  the  Is  state  and  one 
in  the  Ip  state. 

This  seems  to  indicate  that  Ha"  cannot  be  formed  by  collision  of  H  and  H~, 
but  only  by  collision  of  Hg  with  an  electron.  Also,  the  proposed  mechanism  for 
formation  of  H  and  H"  by  dissociation  of  Ha"  may  have  another  step, 

Hs  +  e  —$►  (H2")  — »  (H1B  +  Hap  +  e)  — $>  H  +  H~ 

where  the  Hap  may  recapture  the  electron  of  the  Ha"  or  capture  another  of  the 
electrons  in  the  experimental  system. 


IV  o  Conclusions 

The  potential  curve  and  electronic  wave  function  for  the  resonant  ^g 
state  of  Hg"  have  been  calculated  by  treating  a  one-configuration  trial  function 
as  if  it  were  a  real  variation  function  which  could  he  minimised  to  a  real 
discrete  eigenvalue.  An  average  energy  very  near  to  the  experimentally  ob- 
served energy  was  obtained,  however,  since  this  is  a  state  in  the  continuum 
there  is  no  rigorous  minimum  principle,  and  the  calculated  average  energy  may 
be  above  or  below  the  actual  energy  of  the  system* 

There  is  no  doubt  that  if  the  configuration  determining  parameters  n, 
m,  and  v  were  also  allowed  to  vary,     the  final  result  would  have  been  the 
one  configuration  approximation  to  the  ground  state  energy  of  the  hydrogen 
molecule. 

Excellent  agreement  with  experimental  spectroscopic  constants  has  been 
obtained,  which  substantiates  the  statement  that  this  is  a  resonance  with  a 
lifetime  long  enough  to  be  observable,  and  the  use  of  physical  intuition  to 
pick  a  restricted  variation  function  which  simulates  an  electron  scattering 
off  a  deep,  narrow  potential  well  formed  by  an  excited  state  of  H^o 

The  method  of  finding  the  longest -lived  wave  packet  in  the  continuum 
by  a  restricted  variation  method  coupled  with  good  physical  guesses  may  have 
considerable  application  in  other  problems. 

It  is  the  Intuitive  feeling  of  the  authors  that  all  resonant  states  will 
be  able  to  be  treated  by  this  approximate  method,  because  if  they  have  a  long 
enough  lifetime  to  be  observed,  they  are  probably  weakly  correlated,  in  the 
sense  described  above  and  therefore  amenable  to  approximate  variation  treatment <> 
This  may  not  be  a  necessary  condition,  but  it  is  certainly  a  sufficient  one. 
The  lack  of  rigor  prescribes  great  care  in  any  application* 
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Table  3T 
Wave  Function  and  Potential  Curve 


R  (au)   -E  (au) 


a 


R  (au)   -E  (au) 


a 


1.596  .74336 

1.7196  .74848 

1.9196  .75064 

2.1196  .74804 


1.115    .765 

.445  -.17 

.100  0.0 

1.175   .810 

.485  -.185 

.105  0.0 

1.270    .890 

.550  -.210 

.105  0.0 

1.365  .970 

.600  -,24o 

0I05  d.o 


2.2432  a7hkSS 


3.000    .71513 


7.000    .63377 


10.000   .62929 


15.00 


.61352 


1.420 

.635 
.105 

1.75 
.82 
.11 

3^53 

1.78 

.11 

5.01 

2.51 

.11 

6.24 

3.76 

.10 


1.020 

-.255 
0.0 

1.34 

-.38 

0.0 

3.39 

-1.66 
0.0 

4.96 

-2.50 

0.0 

6.21 
-3-76 
0.0 


a.  Throughout  this  table,  the  electron  parameters  are  listed  in 
the  order  /  0  0  0\  so  that  the  inner  electron  is  first,  and 


/  0  0  0' 

I  001 
\2  0-l, 


the  outer  resonant  electron  is  last. 

b.  This  point  was  not  completely  varied,  but  a  definite  continuation 
of  the  trend  to  Hi8  +  Hip  +  e"  is  exhibited. 


■Cable  II 
Comparison  of  Spectroscopic  Constants 


calculated 

experimental 

E 

-o751^  au 

-,738  au 

«e 

1,9^9  au 

*e 

oJO  ev 

.31  ev 

^eXe 

,042  ev 

*"o01  ev 

Be 

o0039  ev 

Oft 

3  x  1CT4  ev 

D* 

5o36  ev 

Table  III 
Oonparieon  with  H2  ^  Potential  Curve 


R  (au) 

— E  H2C  TTu 

— E  H;»  7L„ 

AS 

1.9 

.71422 

-75057 

.03635 

2.0 

.7l4l8 

.74800 

.03382 

2,25 

.71006 

,74481 

.03475 

h9 

.68445 

.71515 

.03068 

7.0 

.62478 

,63377 

.00899 

10.0 

.62448 

.62929 

.00481 

a.  From  ref.  (12) 


) 
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